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Abstract. We present a search for the hyperfine transi- 
tion of the 57 Fe +23 ion at 3.071 mm in clusters of galaxies 
with the ATNF Mopra telescope. The results are com- 
pared with a realistic estimate of the peak brightness tem- 
perature of the line in a cooling flow cluster A85, using the 
available X-ray data. 

Key words: galaxies: clustering - clusters of galaxies: 
individual Abell 85 - X-rays: galaxies - HII regions: 30 
Dor 



1. Introduction 

It has been suggested by Sunyaev and Churazov (1984) 
that it might be possible to detect the hyperfine transi- 
tion of a Li-like 57 Fe +23 ion line in clusters of galaxies. 
Hyperfine transitions have been very important in astron- 
omy. The first such transitions was predicted for neutral 
hydrogen atoms (the famous 21cm line) by van de Hulst 
in 1945 and detected in 1951. Atomic hyperfine splitting 
is the result of the interaction between the magnetic mo- 
ment of the nucleus and the magnetic moment of the elec- 
tron plus the magnetic field generated by the circulating 
electrons. In the case of the HI 21cm hyperfine transi- 
tion at 1420 MHz, the ground state l 2 Si/2 is split into 2 
sublevels. Deuterium has a similar line predicted at 327 
MHz but there have been no detections to date despite 
many efforts. The spontaneous transition probability (or 
Einstein coefficient) of the HI hyperfine transition is very 
low, W — 2.85 x 10 -15 sec -1 . However, external excita- 
tion mechanisms such as collisions with other H atoms and 
electrons or Lycv radiation increase the transition rates 
enormously. Hyperfine splitting is not present for all atoms 



or ions, since the interaction between the nuclear mag- 
netic moment and the electronic magnetic moment does 
not occur if either of them is zero. We know from quan- 
tum mechanics that the total nuclear angular momentum 
(or nuclear spin) is zero for nuclei with an even number 
of protons and neutrons, and that the total electronic an- 
gular momentum of an ion with completely filled electron 
shells is zero. Thus the commonly available isotope of iron 
56 Fe has no hyperfine structure. However, Sunyaev and 
Churazov (1984; SC) pointed out that the Li-like 57 Fe+ 23 
ion exhibits hyperfine structure with a higher spontaneous 
transition rate (W = 9.4 x 1CP 10 sec -1 ) than that of 
the HI hyperfine transition (W = 2.85 x 10~ 15 sec -1 ), 
which compensates the low abundance of the 57 Fe isotope 
( 57 Fe/ 56 Fe=2.3%; Vokening 1989). The wavelength of the 
57p e +23 n yp er fi ne transition has been calculated theoreti- 
cally to be 0.3071(3)cm (Shabaeva & Shabaev 1992). The 
FWHM of the line is 40 MHz (or 130 kms -1 ), assuming 
a cluster gas temperature of 2 x 10 7 K and that the line 
width is due to thermal Doppler broadening alone. 

Clusters of galaxies are known to possess hot gas 
of temperature 10 7 — 10 8 K. This intracluster medium is 
found to be metal rich with typically half solar abundance 
(Sarazin 1988). The early X-ray observations detected the 
prominent "7 Kev Fe line" which is a blend of Fe +24 and 
Fe +25 Ka and Kf3 lines. The detection of this line not 
only provided the concrete proof that the primary X-ray 
emission mechanism for clusters is thermal bremstralung, 
but also showed that the intracluster gas has nearly solar 
abundances and that at least some of the gas must have 
been ejected from galaxies. Other Fe lines at lower energy 
such as Fe +23 L lines have been detected at ~ 1 keV in 
clusters. Thus we know that Fe +23 is abundant in cluster 
environments. SC had estimated the 57 Fe +23 line strength 
in the centre of the Perseus cluster to be ~ 1.5 mK, as- 
suming a central temperature of 1.4 keV inside a 10-20 
kpc radius with an X-ray emission integral of ~ 4.8 x 10 66 



Send offprint requests to: H. Liang, h.liang@bristol.ac.uk 



cm 



-3 



2 



H. Liang et al.: Observations of the 57 Fe +23 Hyperfine Transition in Clusters of Galaxies 



The detection of the 57 Fe +23 line will have significant 
astrophysical implications on our understanding of the 
properties of the galaxy clusters. Firstly, it will provide 
us with a direct measure of the dynamics of the ICM, 
such as turbulence and accretion flows. There have been 
great debates on the formation of tailed radio sources in 
clusters of galaxies (e.g. Burns et al. 1986; O'Donoghue 
et al. 1993; Roettiger et al. 1993), with some suggesting 
that the sources are bent by their motions through the 
ICM and others suggesting that the bending is the result 
of an intracluster wind. This debate would benefit greatly 
from a direct measure of the intracluster turbulence. Sec- 
ondly, mm observations of the 57 Fe +23 line will comple- 
ment the X-ray measurements of the distribution of iron 
in clusters and measure the 57 Fe/ 56 Fe ratio which would 
determine the nucleosynthesis history of the intracluster 
gas. Finally, the mm line would provide a direct and high 
precision measure of the cluster redshift since the spec- 
tral resolution at mm-wavelength is much higher than in 
X-rays. 

2. Observations 

Our search for the A 3.1-mm 57 Fe +23 line was carried out 
in 1996 September using the Mopra telescope of the Aus- 
tralia Telescope National Facility. The antenna is 22m in 
diameter of which the inner 15m is illuminated by the 3- 
mm feed, giving a beam width (FWHM) of 45" and gain 
of 40 Jy K . The receiver is a cooled SiS mixer with 
receiver temperature of 80 K which was checked with hot- 
cold load measurements twice during the seven day obser- 
vation period. Pointing was tested by using SiO masers, 
the largest pointing error measured was 10" . Atmospheric 
opacity was checked by skydips, the zenith opacity varied 
between 0.06 and 0.12 for all but a few hours during the 
run. 

The spectrometer is an autocorrelator with two sec- 
tions, each covering 256 MHz with 256 channels (1 MHz 
per channel). We overlapped these two by 125 MHz so 
as to cover a total bandwidth of 384 MHz. Roughly 50 
MHz on each side are compromised by the IF filter shape, 
leaving a total bandwidth with consistent sensitivity of 
about 285 MHz which translates to 920 km s _1 at a wave- 
length of 3.2 mm (the redshifted wavelength). Since the 
line width is not well known, we do not remove baselines 
from the spectra, but simply edit the records on the basis 
of their raw rms noise. 

Each record consists of a pair of two minute on and off 
source spectra from which we take the difference normal- 
ized by the bandpass shape. The reference position is offset 
by six minutes of right ascension. The expected spectral 
rms per record is about 20 mK, which is what we find in 
almost all cases. A few records (<~ 10%) show significantly 
higher noise due to baseline instabilities in the IF system. 
They are removed from the average. With this editing, 
the rms noise in the average decreases as integration time 



s/Unt up to our maximum integration time of 10.5 hours. 
Because the spectra have not been baseline fitted, there 
remains some structure in the average spectra, but it is 
quite broad band (typically 100 MHz wide or greater). 
This baseline structure limits our ability to detect broad 
lines. 

The clusters which we are searching for the 57 Fe line 
are at redshifts of ~ 0.06. Assuming rest frequency of 
97.611 ± 0.03 GHz (A = 3.0713 mm) the cluster center 
frequencies come out in the range 92 - 92.5 GHz. More 
nearby clusters would not necessarily give stronger sig- 
nals, since the beam size of 45" translates to only 35 h _1 
kpc for a cluster at a redshift of 16,000 km s _1 . Thus the 
emission region is probably considerably larger than the 
beam. This is the reason for taking the reference position 
so far from the source. 
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Fig. 1. Spectrum of the H41a line from 30 Dor. The antenna 
temperature is plotted against the helio-centric velocity for the 
data from both spectrometer channels separately and com- 
bined. The vertical separations of 30 mK are for clarity in 
display and arbitary. 

As tests of the sensitivity and stability of our system 
and observing strategy, we have observed several calibra- 
tion lines. The first test is simply to check the pointing 
using SiO masers at 86 GHz. We did these pointing checks 
every few hours for the first several days of the run ob- 
serving five positions offset by a half beamwidth from the 
tabulated source position. We always detected the masers 
within 10" of their nominal positions and with their nomi- 
nal line strengths. Since both the beam size and the cluster 
sizes are much larger than this we have not introduced any 
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pointing offsets. This test is limited in that the 86 GHz 
SiO frequency requires a different receiver tuning than the 
92 GHz 57 Fe line, so we checked the receiver performance 
by observing test lines at 92 GHz. The best candidates are 
13 CS line from Orion (SiO position) and M17SW, whose 
frequency is 92.494 GHz and the H41a recombination line 
at 91.955 GHz which we observed in the 30 Dor region of 
the Large Magellanic Cloud. The 13 CS line is so strong 
that it shows up in a single on-off pair, but the H41a line 
is weak enough to provide a good test of a moderately 
long integration. Our spectrum of this test line is shown 
in Figure 1 ; the position chosen was the 20-cm continuum 
peak at 05:38:46.675, -69:04:58.6 (J2000). The total inte- 
gration time is 112 minutes and the noise level is 2.5 mK, 
2.0 mK would be predicted by the radiometer equation, 
so the system is evidently performing well for integration 
times of an hour or so. 




Table 1. Summary of Observations 
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mK 10 2 km s _1 Fig. 2. Spectrum for the cluster A85. Antenna temperature 
is plotted against the helio-centric velocity for the 2 spectrom- 
eter channels separately. The vertical separation of 6 mK is 
arbitary, for clarity in display. 
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Table 1 lists the clusters which we have searched for 
57 Fe emission. Columns 1-3 give the cluster name and po- 
sition, column 4 gives the integration time on-source, col- 
umn 5 gives the detection threshold, and column 6 gives 
the useful velocity range. The detection threshold is a 
strong function of the width of the line because we have 
subtracted no baselines from our spectra. The threshold 
given on Table 1 applies to lines narrower than about 150 
km s _1 in full- width. For broader lines the residual base- 
lines in the receiver system begin to degrade our detection 
threshold. This effect is illustrated in Figure 2, which is 
the pair of spectra for the cluster A85 for the two spec- 
trometer channels. The vertical offset between the two is 
arbitrary, here it has been set at 6 mK for convenience. 
One channel is dominated by a linear baseline, while the 
other shows primarily a quadratic baseline shape. If we 
were to remove first and second order baselines and align 
the two channels we would obtain a spectrum with no 
features above 2 mK or less (depending on the smooth- 
ing). On the other hand, if cluster lines are as broad as 
several hundred km s _1 then this baseline removal could 
remove a real feature, and so our threshold would have 
to be raised; we cannot rule out the lump between 16200 
and 16500 km/s being part of a real feature of strength 
3 mK and lincwidth ~300 km s _1 . For lincwidths of 200 
km/s or greater we must set the threshold at 5 mK in 
this case. For the other clusters on table 1 there are simi- 



lar baseline irregularities that set the maximum lincwidth 
corresponding to the threshold at about 150 km/s. 

3. Discussions and Conclusion 

Following the formalism given in Sunyaev and Churazov 
(SC;1984), we can calculate the expected line intensity for 
a cluster given the cluster gas parameters derived from 
X-ray observations. 

The peak brightness temperature of the 57 Fe +23 hy- 
perfine transition should be 



AT b = ^ - 1 a g "fr 
8tt 3 / 2 ^ k /S.vd gi 



(1) 



where g u ,gi are the statistical weights of the hyperfine sub- 
levels, Ais£> is the Doppler line width (i.e. the half-width 
at e _1 of the maximum), a is the cluster isotope abun- 
dance relative to hydrogen and the function $ is given by 



/ 



WS(T)W e (av) 



W{1 + N) + \n e {av) + ^WN + ^\n e {av) 



dl{2) 



where W is the spontaneous transition probability, S(T) 
is the fraction of ions of a given species at temperature T, 
N is the occupation number of the microwave background 
at the wavelength of transition, (av) is the excitation co- 
efficient by resonance scattering which is the dominant 
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excitation process (SC). Note that Eq. (||) differs slightly 
from equation (2) in SC (see Appendix). 

In the case of clusters of galaxies where WN ^> 



^n\{av) (or n e <C I cm 3 ), Eq. (2) is reduced to 



4> 



1+N+^N e 

91 



(3) 



Assuming a constant temperature, we have from Eq. (1) 
and Eq. (3) 



ATh oc / n dl oc S x 



(4) 



where S x is the X-ray surface brightness. If we assume 
solar abundance for clusters, where solar Fe abundance 
was taken to be 8.8 x 10 -7 (SC), and isothermal cluster gas 
temperature within the Mopra beam to be at its optimum 
value of 2 x 10 7 K where S(T) has its peak value of ~ 
0.3 (Jacobs et al. 1977), we can deduce a simple relation 
between the peak brightness temperature of the 57 Fe line 
and the X-ray surface brightness: 



AT h ~ loos',. 



(5) 



where AT), is the peak brightness temperature of the line 
in units of mK before convolution with the beam of the 
radio telescope, and S x is the X-ray surface brightness 
in units of ergss -1 cm~ 2 ster _1 in the ROSAT band of 
[0.1,2.4] keV after corrections for the telescope response 
and Galactic absorption. In the case of an average clus- 
ter of galaxies, where n e ~ 10~ 3 — 10~ 2 cm~ 3 , the peak 
brightness temperature would be ~ 0.001 mK. However 
in the centre of cooling flow clusters of galaxies where the 
central density could be as high as n e ~ 0.1 — 0.2 cm~ 3 , 
the peak brightness temperature could only reach ~ 0.1 
mK in the most extreme cases, contrary to the estimate 
of ~ 1.5 mK given by SC for the Perseus cluster. 

The cluster A85 is considered a classic example of a 
cooling flow cluster with high X-ray luminosity and thus 
one of the best candidates for this experiment. To obtain 
a realistic prediction for the 57 Fe line brightness for A85, 
we examine the X-ray data available for the cluster. The 
ROSAT PSPC data for A85 has been analysed by Pislar et 
al. (1996), who found an X-ray luminosity of 9.3±0.2x 10 44 
ergs s~ 4 within a radius of 1.4 Mpc in the ROSAT band of 
[0.1-2.4]kcV. The X-ray surface brightness was fitted with 
a [3 model given by, 



oc [1 + (r/r c ) 2 r 3 ^ +1 / 2 
which corresponds to 
n e (r) = n Q [l + (r/r c ) 2 ]- 3 ^ 2 



(G) 



(7) 



for the electron density distribution, where uq = 
0.01cm- 3 , (3 = 0.438 and r c = 32kpc (25"). However, 
the surface brightness profile obtained with the ROSAT 



HRI (Prestwich et al. 1995) which has a resolution of 5" 
(c.f. 20" the PSPC) clearly shows that there is an excess 
above the j3 model within the central 20" radius. Taking 
the X-ray surface brightness from HRI, we obtain from the 
above relation a peak brightness temperature of 0.003 mK 
for the 57 Fe line. This predicted line brightness can be in- 
creased, if the abundance of 57 Fe is significantly different 
from the terrestrial abundance and/or the distribution of 
Fe increases towards the centre. However, such effects are 
unlikely to increase the detectability of the line by more 
than a factor of 10. 

This analysis shows that the expected brightness tem- 
perature of the 57 Fe line to be <C 1 mK in cluster cores, 
even for an archetype cooling flow cluster such as Abell 
85. We have shown that the MOPRA telescope can reach 
a sensitivity level comparable to the original prediction 
given by SC of the line intensity but is certainly inade- 
quate given the more realistic estimates shown above. 
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Appendix: derivation of Eq. (1) & (2) 

The number of transitions upwards equals to the number 
of transitions downwards: 



n u W + n u B u iB u + n u c u i = niB lu B v + 



(i) 



where n Ul rii are the number in the upper and lower 
states respectively, W and B are the Einstein A and B 
coefficients relating to spontaneous and stimulated emis- 
sion/absorption, B v is the intensity of the microwave back- 
ground and c is the collisional excitation/de-excitation 
rate. The W,B and c coefficients are related as follows 
(Rybiki and Lightman 1979): 

Bin = 9 —B uh B ul B v = WN, c lu = ^c ul eM~hv/kT)(2) 
9i 9i 

where N = (expQiv/kT) — 1) _1 is the occupation number 
of the microwave background. As was discussed in SC, the 
dominant excitation mechanism in this case is resonance 
scattering of an electron by the ion and the rate of excita- 
tion is given by c u i — \n e (&v) , where (av) is the average of 
the product of the resonance scattering cross section and 
the relative velocity between the electron and the ion. 
The peak intensity of the line is given by 

h = {N U W + N u B ul B v - NiBi u B v )^-<l>(yo) (3) 

47T 

where N u ,Ni are the total number of ions in the upper 
and lower levels per unit area along the line of sight re- 
spectively, W is the spontaneous emission coefficient or 
Einstein A-cocfficicnt, is the central frequency of the 
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line and <j> is the line profile. Assuming thermal dopplcr 
broadening alone, we have 

<t>{y) = -;Jr e -(--"o) 2 /A^ f4) 



where Ai/jj = y jj^vq is the doppler width which is re- 
lated to the FWHM Avl = 2\/ln2Az/£). In our case where 
hv <C kT, we have the peak brightness temperature of the 
line given by AT;, = ^^loi^o)- By substituting Eq. 1 into 
the expression for Iq in Eq. (3), we have 

AT b = W |^ M - Nu \n e{ * V )) (5) 

where 

N u = f a5(T) Hu —n e dl (6) 
J ni+n u n e 

with ^ ~ 1 and a5(T) giving the abundance of 57 Fe +23 
in the cluster, and Ni — (ni/n u )N u . From Eq. (I) and (2), 
we have 

nu _ g u WN+ \n e {av) 

ni ~ gi W{N + l) + \n e {av) 1 ' 

Thus the peak brightness temperature is given by 

87r ci /' 2 ^ k Avd gi 

where 

$ = / v ;2 eV : dl(9) 

J W{l + N) + \n e {av) + ^WN + ^\n e {av) 

The term ^■■^n e (av) in the denominator is absent in equa- 
tion 2 of SC, however in cases where WN ^> ^nl(av) this 
term is insignificant. 
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